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abstract 

\ We present the results of a deep optical 2.6 square degree survey with near infrared 

follow-up measurements of the intermediate- aged Praesepe open cluster. The survey 
is complete to Z c =21.3, 2=20.5, corresponding to 0.06 M Q assuming a cluster 
age of 0.5 Gyrs. Using 3-5 pass-bands to constrain cluster membership, we identify 
32 new low mass cluster members, at least 4 of which are likely to be substellar. We 
use the low mass census to trace the region where the sequence moves away from the 

■ NEXTGEN towards the Dusty regime at T e ff = 2200K. In doing so, we identify four 
' unresolved binaries, yielding a substellar binary fraction of ~ 30 percent. The binary 
| fractions appear to decrease below 0.1 M©, in contrast to the rising fractions found 

in the Pleiades. We also identify a paucity of late M dwarfs, thought to be due to a 
steepening in the mass-luminosity relation at these spectral types, and compare the 
properties of this gap in the sequence to those observed in younger clusters. We note 
an overdensity of faint sources in the region of the so-called subcluster (possibly an 
older smaller cluster within Praesepe), and subsequently derive the luminosity and 
mass functions for the main Praesepe cluster, revealing a turn-over near the substellar 
boundary. We conclude by presenting astrometric measurements for low mass Praesepe 
candidates from the literature, and rule out as a likely foreground dwarf RPrl, hitherto 
thought to be a substellar member. 
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1 INTRODUCTION 

There has been much work done surveying young open clus- 
ters for ve ry low mass stars (VLMS) and brow n dwarfs 
(BDs) (e.g. lBouvier et alll998LlDobbie et"alll2002l) . The fo- 
cus has generally been to identify BDs and measure clus- 
ter mass functions (MFs) in a variety of environments 
to establish if the initial mass function (IMF) is univer- 
sal or not - the answer to which has important impli- 
cations for our understanding of star and BD formation. 
The study of older clusters is also important if we are 
to understand how the dynamical evolution of the clus- 
ter affects the shape of their MF - over time VLMS 
and BDs are expected to be preferentially eje cted (eg. 
Ide la Fuente Marcos fc de la Fuente Marcos! feoOCft . For the 
Hvades. iDobbie et alJ l|2002f) searched for low-mass (0.1- 
0.06Mq) members in a 10.5 square degree survey, but found 
only one (previously known) stellar member. They estimate 
that 4-5 Hyads should have been found in this mass range 
if the Hyades and Pleiades MFs are identical, and conclude 



that dynamical evaporation of low-mass Hyads is the most 
likely explanation for the deficit. However, studies of other 
clusters with similar ages are clearly needed if we are to 
properly address this question. 

Late cluster members are also ideal for testing ultra-cool 
atmos phere models of late M and L spectral types (T e // ^ 
2500K lLeggett et aljfeoOlT) . At such low T e// , atmospheric 
dust grains condense out of the gas phase, and strongly af- 
fect both colour and spectral properties. However, dust grain 
properties depend not only on temperature but on gravity 
(g)and metallicity ([M/H]) which cannot be measured with 
confidence for late M and L field dwarfs. But with well con- 
strained age (and hence radii and masses inferred from evo- 
lutionary models) late M and L cluster members' g, T e // and 
[M/H] can be known. Indeed, spectroscopic studies of such 
objects with a range of different (but known) ages would 
allow us to empirically isolate the Teff, g, and [M/H] de- 
pendence of both narrow and broadband spectral features. 

The Praesepe cluster has an age of (0.9±0.5) Gyrs, 
lies at a distance of ~ 170 pc, has a near so- 
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lar metallicitv llHambly et all Il995ai and zero redden- 
ing llCrawford fc Barn^a^^69jT. Large scale proper motion 
studies of Praesepe include Irlamblv et alJ l|1995bll (here- 
after HSHJ) coveri ng 19 sq. degs down to 0.1M©, and 
lAdams et all (2002) covering the whole cluster to similar 
depth. A number of small er but deeper Praes epe surveys 
have also been carried out. iPiniield et al] dl997l) (P97 here- 
after) covered ~ 1 sq. deg down to 7 C = 21. Ipinfield et ail 
l20Q0j) (POO hereafter) covered ~ 5 sq. degs down to 7~19.5. 
Mag azzu et alJ lll998li covered 800 sq. arcminutes down to 
/=21.2, and discovered one candidate member of M9 spec- 
tral type (Roque Praesepel; RPrl hereafter). If confirmed as 
a Praesepe member, its estimated mass would be M=0.063- 
O.O84M0 making it a possible BD. NIR cha racterisation of 
the fa int c andidates has been p resented bv iHodgkin et ail 
lll999l) and lPinfield et all d2003li (P03 hereafter). 

In this paper, we present the results of a deep 2.6 square 
degree optical survey of the cluster, and use near-infrared 
follow-up measurements to refine photometric membership 
status. We then combine our candidates with others from the 
literature, and identify likely unresolved binaries. We also re- 
cover the "M dwarf gap" in the cluster - a dearth of M7-8 
dwarfs previously n oted in the Pleiades and other clusters 
iDobbie et aljl2002l : P03). We derive the luminosity func- 
tion, use state-of-the-art evolutionary models to estimate 
candidate masses and determine the cluster mass function. 
We also derive proper motions for several previously identi- 
fied cluster candidate members. Finally, we discuss planned 
future work. 



2 WFC IZ SURVEY 

2.1 Observations 

Deep / and Z band images of Praesepe were obtained with 
the INT/WFC during the nights 2001/12/24-28. The WFC 
instrument consists of 4 thinned 2048 x 4196 pixel CCDs 
(0.333 arcseconds/pixel), and operates at the prime focus of 
the INT, covering a 0.29 sq. deg field of view. 

The conditions for ~ 35 percent of the time were pho- 
tometric with modest humidity and seeing of ~ 1 arcsec. 
During the run, we observed 9 pointings through the I rgo 
and Z rgo filters with exposure times of 20 minutes per band. 
In total, 2.6 square degrees of the cluster were surveyed. Our 
coverage is shown in Figure 

2.2 IZ reduction 

Data reduction was performed using standard IRAF 1 rou- 
tines. The images were bias subtracted, and non-linearity 
accounted for prior to flat-fielding. Fringe maps, which were 
constructed in both bands by median filtering images ac- 
quired throughout the entire run, were used to remove the 
effects of the interference between night sky lines in the CCD 
substrate. 

Sources were identified in the /-band images using 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of of Universities 
for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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Figure 1. INT survey sky coverage. The 9 WFC fields are shown. 
2 circles of radii 0.5 and 1.5 deg centred on the cluster centre are 
also shown. Point sources from the USNO catalogue are shown as 
filled circles, with size indicating apparent magnitude. 

IRAF'S DAOFIND at a detection level of 3a, and matched 
to Z- counterparts using a 4 pixels search radius of the I- 
band source position (taking into account any systematic 
shifts between the /- and Z- images) . Aperture photometry 
was measured for each source using the PHOT routine with 
a 7 pixel aperture size and a sky value computed from within 
a circular annulus. 

To calibrate the photometric data, standard stars drawn 
from lLandoltl 1119921) were observed, and zero points and ex- 
tinction coefficients determined. In order to calibrate the 
7 images obtained in non-photometric conditions, we re- 
observed these fields with shorter (45s) integrations in the 
I- and Z- bands on the night of 23rd February 2002 dur- 
ing photometric conditions. Landolt standard observations 
were used to calibrate the magnitudes of ~ 50 of the brighter 
stars in each of these shorter integrations, and of these, the 
~ 40 stars that were not saturated in our longer integrations 
were used as secondary standards to determine zero points 
for these longer exposure images. 

2.3 WFC Results 

Figure |21 plots the frequency of the total number of sources 
detected per 0.1 magnitude interval in the I- and Z-bands. 
The dotted lines are straight line fits to the I- and Z- band 
histograms in the magnitude ranges /=17-20.5 and zT=17- 
19.5. Survey completeness limits for the / and Z bands were 
taken to be the magnitudes at which the binned star counts 
drop below 90 percent of the expected frequency, and are 
shown as vertical dashed lines. In this way, the survey was 
estimated to be 90 percent complete to / =21.3, Z—20.5. 

Our WFC survey resulted in a catalogue of ~ 15000 ob- 
jects. The resultant I,I-Z colour-magnitude diagram (CMD) 
is shown in Figure El In order to separate candidate 
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Figure 2. / and Z magnitude histograms of source counts in 
the Praesepe WFC survey. I- (bold) and Z- band data are shown 
by solid lines. Dashed lines indicate the 90 percent completeness 
limits. 



cluster members from background sources, it was neces- 
sary to locate the cluster sequence. To do this, we used 
0.5 Gyr DUSTY JChabrier et al.1 ExioT) and NEXTGEN 
fearaffe et alj|i"99£ ) model isochrones (dashed and dotted 
lines respectively), transformed onto our CMD. Predicted 
I-Z colours were c alculated using the relations given in 
iDobbie et alJ l|2002l) . where we assume a distan ce modu- 
lus of (m-M)o = 6.16 jPinsonneault et al]ll99ct) and zero 
extinction. To further highlight the location of the cluster 
sequence, we also identified HSHJ proper motion members 
and optical/NIR candidates from P03 that were recovered in 
our survey. We choose not to include objects that fall within 
1 degree of the centre of a suspected sub-cluster, discussed 
in &\ These are highlighted in the figure as plus signs and 
triangles respectively. Further, we overplot candidates from 
P03 that are shown to be astrometric cluster members in 
section §10\ and are shown as filled stars. 

Also plotted in Figure [3] are the expected magnitudes 
of several M field dwarfs as they would appear at the clus- 
ter distance. Praesepe members are sufficiently old to have 
finished their contraction stage, implying that cluster UCDs 
will have the same luminosity as field stars o f the same T e f f . 
The m agnitudes for these are taken from iGossburn et alJ 
( 1998), and are shown as squares in the CMD. 

Using the photometric and astrometric members, com- 
bined with the offset field star positions and the synthesized 
model tracks as a guide on the CMD, we have defined a 
cut-off line to separate background sources from potential 
cluster members. By taking into account photometric un- 
certainties (indicated on the right hand side of the figure), 
we have positioned our cut-off line such that it is at least 
1-sigma (in I-Z uncertainty) blueward of these photometric 
guides. Note that for I> 19, the offset field stars are the clos- 
est objects to the cut-off line, and the proper motion clus- 
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Figure 3. WFC IZ colour magnitude diagram. The solid line 
indicates the redness cut-off for Praesepe members. The dashed 
and dotted lines show the theoretical isochrones for the 0.5 Gyr 
DUSTY and NEXTGEN models respectively. Plus signs and tri- 
angles are cluster LMS from HSHJ and P03 surveys respectively. 
Proper motion members identified in this work ar e shown as filled 
stars . The expected positions of field dwarfs (from lCossburn et alJ 
Il998l) are shown as squares. 

ter members are rather redder than the cut-off. This is as 
expected, since Praesepe low-mass members (Teff<2500K) 
will not be fully contracted and we expect lower Teffs as a 
function of magnitude (e.g. 100K lower for Mi 14 at 0.5Gyr 
age), and thus slig htly larger I-Z colours (e.g. 0.1 larger; 
Dob bie et all 12002). The use of the offset field stars thus 
provides a conservative way to define our cluster region on 
the CMD. We are thus confident that our selection will not 
miss any genuine cluster members. We further constrained 
the selection of candidates using a bright limit (brighter 
than which 2MASS is sensitive to clu ster members which 
will have been previously identified bv lAdams et al.l ll2002l) 
and HSHJ), and a faint completeness limit (both indicated 
on the CMD). All stellar-like objects redward of the cut-off 
line, and between the bright and faint limits were selected as 
IZ candidates, and are shown as filled circles in the figure. 
We identified 320 cluster candidates in this way from J c = 
17.5-21.3. Visual inspection of each revealed 95 to be spuri- 
ous (mostly diffraction spikes and cosmic ray hits), leaving 
225 genuine candidate members. 



3 NIR FOLLOW UP 

Near-infrared follow-up measurements of the 225 WFC can- 
didate members were made using the Fast Track Imager 
(UFTI) instrument on the United Kingdom Infrared Tele- 
scope (UKIRT). UFTI was windowed (512x512 pixels) for 
faster readout providing a 46.5 arcsecond field of view, with 
each integration consisting of a 5-point dither pattern. K- 
band photometry of the 120 brightest candidates was mea- 
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sured during flexibly scheduled queue observing in photo- 
metric conditions from the 19th to 24th of February 2003. 
K-band photometry of the remaining 105 fainter candidates, 
as well as longer J-, H- and K-band integrations were ob- 
tained by us in photometric conditions during our UKIRT 
observing run on 7-10 March 2003. 

Our observing strategy was to initially measure ri- 
band photometry (S/N~20) for all candidates, with expo- 
sure times ranging between 6s to 60s depending upon source 
brightness. At this stage we considered a target as a poten- 
tial cluster member if it lay near to the NEXTGEN isochrone 
for Jf<15.5, had I-K>3.3 for #=15.5-16.5, or had I-K>3.8 
for K>16.5. We then obtained longer integration J, H and 
K measurements of 12 of the faintest potential cluster mem- 
bers, where we used exposure times of 24s-120s, aiming for 
a S/N of ~ 50. The optical and NIR photometry of all WFC 
sources flagged as potential candidate members is presented 
in Table |A"T1 

The images were de-biased, dark subtracted, flat fielded 
and finally combined into mosaics using t he UKIRT software 
pipeline ORACDR ferideer et aUl200Ch . The STARLINK 
package GAIA was then used to measure the photometry, us- 
ing apertures matched to the seeing (0.6-1 arcseconds). In- 
dividual aperture corrections were calculated using UKIRT 
faint standards measured either side of the candidate, and 
instrumental magnitudes were then transformed onto the 
MKO system using airmass curves derived from these stan- 
dards. 



4 ADDITIONAL PRAESEPE CANDIDATES 

As well as the new WFC candidates, we also consider Prae- 
sepe candidate members from several other surveys. Proper 
motion candidates from HSHJ define the bright main se- 
quence down to 7 C ~ 17.5. We have transformed the HSHJ 
photographi c In magnitudes onto the Cousins system us- 
ing |Bessei3 (Il986f) . Sources from P97 and POO (which we 
will refer to as Riz and Iz candidates respectively) were also 
considered in the present analysis. NIR photometry of the 
Riz and Iz candidates was taken from P03, and is on the 
MK O system. We als o considered the RPrl candidate from 
iMaeazzu et~afl Jl998t) (the M98 survey hereafter). The NIR 
photometry of this object has been trans formed onto the 
MKO system using the transformations of lHawarden et all 
feOOll) . 

In order to obtain NIR measurements of as many of 
the HSHJ stars as possible, we cross matched their posi- 
tions with the 2MASS All Sky catalog. 2MASS photometry 
typically provides a SNR~10 for J = 15.8, H — 15.1, and 
Ks — 14.3. We identified NIR counterparts within a search 
radius of 3 arcseconds from each source. Of the 459 HSHJ 
candidates that have astrometric membership probabilities 
> 70 percent (from HSHJ), 447 had a corresponding 2MASS 
counterpart. 3 of these were flagged by 2MASS as having po- 
tentially contaminated photometry. These, and the 12 un- 
matched HSHJ candidates, are noted in Table \K2\ We trans- 
formed the 2MASS photo metr y onto the MKO system using 
equations in lCarpenterl i200ll) . 



5 PHOTOMETRIC MEMBERSHIP AND 
BINARITY 

5.1 Membership from the IK CMD 

Figure 0] shows the K,I-K CMD for Praesepe. Iz and Riz 
candidates from P97 and P03 are shown as triangles and 
squares respectively. The new WFC candidates are shown 
as diamonds. Large diamonds represent WFC candidates 
that were observed using longer K, J and H integrations. 
All candidates with photometry which we deem inconsistent 
with cluster membership (see below) are displayed as empty 
symbols. RPrl is shown as an asterisk. 

Over-plotted on the CMD are 0.5 Gyr NEXTGEN 
(solid line) and Dusty (dotted line) model isochrones from 
the Lyon Group tearaffe et al]|l998l Iciiabrier et alJl2000l) . 
where we have assumed a distance modulus of (m-M) = 
6.16 and zero reddening. 

The large number of HSHJ and Iz sources define the 
sequence well down to 7\~15.5, and it is clear from the fig- 
ure that the NEXTGEN isochrone fits the cluster sequence 
fairly well down to this magnitude (out to I-K~3.3; as was 
noted by P03). We therefore adjudged the membership sta- 
tus of all candidates with I-K ^ 3.3 based on their proximity 
to the NEXTGEN isochrone. Candidates which appeared to 
be too blue, de spite making allowa nce for the cluster depth 
(+/-0.15 mags; Holl and et aife oOOl and the photometric un- 
certainties, were rejected. 

All 444 HSHJ objects with 2MASS photometry have 
I-K < 3.3. Using the NEXTGEN isochrone as a test for 
cluster membership, we adjudge 18 of these to be non- 
members. Furthermore, three HSHJ candidates appear sig- 
nificantly redder than the isochrone, with photometry in- 
dicative of foreground field dwarfs. However, we caution that 
these sources are all found on the same Schmidt plate in 
HSHJ suggesting that their redness could be due to calibra- 
tion systematic errors, so we choose not to reject them as 
candidates. In Table IX2l we list these 3 HSHJ candidates, 
the photometric non-members and the HSHJ candidates for 
which no 2MASS counterpart was found or a 2MASS con- 
tamination flag was logged. All other HSHJ candidates have 
2MASS photometry consistent with cluster membership. 

For I-K > 3.3 (T e// < 2500K) the cluster sequence 
is not obviously defined by the candidates themselves. Also 
it is not immediately apparent if either the NEXTGEN or 
DUSTY isochrones are appropriate. This is the T e ff range 
when atmospheri c dust grains are exp ected to condense out 
of the gas phase Jjones fc Tsuiilll997Tl . so we might expect 
the NEXTGEN isochrone to be inadequate. However, the 
DUSTY isochrone is clearly too blue when one moves fainter 
than #=15.5, and this isochrone is thus not a viable alter- 
native to guide candidate member selection. Consequently, 
we chose not to rely on the models to define our Praesepe 
sequence in this range. 

Instead, we defined a cluster sequence region in which 
to identify candidates as likely members. In order to do this 
we began by defining the top of the equal mass unresolved 
binary sequence in the K,I-K CMD, which we assume con- 
sists of the brightest candidates in this colour range - these 
sources then define the bright limit of our cluster sequence 
region. The single star sequence will be 0.75 magnitudes 
fainter than the equal mass binary sequence. We allow for 
the broadening of both the single and binary sequences by 
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Figure 4. The Praesepe K,I C CMD. Filled symbols represent cluster candidates (binary candidates are circled), and empty symbols 
are photometrically rejected candidates. Circles are HSHJ candidates, squares are Iz candidates and triangles are Riz candidates. For 
the current survey, diamonds are K band only WFC candidates, and large diamonds are JHK WFC candidates. The solid and dotted 
isochrones are 0.5 Gyr NEXTGEN and Dusty models respectively. RPrl is shown as an asterisk. The greyed area and dashed line within 
defines the cluster sequence region, and is discussed in the text. 



cluster depth effects and assume photometric uncertianties 
of +/-0.2 magnitudes are applicable to both single and bi- 
nary members. Therefore, as illustrated by the shaded zone 
in Figure [I] our bright and faint selection criteria are sepa- 
rated by 1.45 magnitudes. We adjudge candidates with I-K 
> 3.3 lying outwith this region as non-members. 

To summarise, of the 225 WFC candidates with if-band 
measurements, 90 look like non members, and 36 have I-K 
colours consistent with cluster membership (84 percent con- 
tamination amongst the candidates from the I, I-Z CMD). 
Of the 12 WFC candidates for which we measured long in- 
tegration J, H, and K magnitudes, 3 are now flagged as non 
members, because although red (and initially considered in- 
teresting) they lie below our cluster sequence region in the 
CMD. Our IK membership criteria are given in Table \KT\ 
and all IK members (including those from the other surveys 
considered) are listed in Table lA"3l 



5.2 Identifying Unresolved Binarity 

Possible unresolved cluster binaries were identified in the /, 
I-K CMD by dividing our cluster sequence region in two, 



using a division that is mid-way between the region's bright 
and faint limits. This division is shown in Figure 4 as a 
dashed line. Candidate cluster members above this line will 
lie closer to the equal mass binary sequence than to the sin- 
gle star sequence, and have been duly flagged as binary can- 
didates (these sources are circled in Figure 4). Unresolved 
binarity amongst the Riz and Iz candidates (as well as the 
HSHJ sample) has been previously addressed by P03. We 
have nothing significant to add to their analysis of candi- 
date members with I-K<3.3. However, our current sample 
of redder candidates is larger and goes fainter than P03. 
Therefore we have re-assessed the unresolved binarity of any 
Riz and Iz candidates with I-K^3.3 (as well as the WFC 
candidates) using our current approach. Unresolved proba- 
ble binaries are flagged in the last column of Table A3. 



5.3 The JK CMD 

The K,J-K CMD is shown in Figure All faint candidates 
with accurate J-band measurements are plotted. We only 
include HSHJ sources with K- uncertainties <0.06 magni- 
tudes, and we do not show Iz candidates with photonic- 
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Figure 5. The Praesepe K, J-K CMD. Symbols are the same 
as in figure |3] Th e M- L spectral boundar y, estimated from 
ILeggett et alj ll998l) and iKnapp et alJ l2004) , is indicated by a 
dashed line at J-K=1.2. WFC candidates discussed in the text 
are labelled. 



try from 2MASS due to the large photometric uncertainties. 
NEXTGEN and DUSTY model isochrones for 0.5Gyr age 
are shown as solid and dotted line respectively. It can be 
seen that the NEXTGEN isochrone agrees well with the 
cluster members to K ~ 15.5 (as with the IK CMD). How- 
ever, for K>15.5 the candidates move significantly red-ward 
of the NEXTGEN isochrone. In the main, these fainter can- 
didates lie slightly above the DUSTY isochrone. However, 
the NIR colour trend shows the rapid increase in J-K that 
the DUSTY models predict, although it occurs at a slightly 
brighter magnitude. 

Three candidates with J-band photometry were as- 
signed as non-members in §5 (two of which appear on the 
plot as open symbols). WFC30 and WFC22 lie slightly below 
the dusty isochrone but have J-K colours that are compa- 
rable with the majority of faint candidate members. This is 
consistent with these two sources being background stars. 
WFC57 has a very red colour ( J-K— 1.91) and lies redward 
of the plotting area. Although this is suggestive of a late L 
dwarf, the I-K colour of 3.88 is not red enough to be a late L 
dwarf (cf. I-K =4-5). Consequently we believe this source is 
a red galaxy. Izl is also highlighted in Figure^ and appears 
to look like a bright foreground field dwarf, but the possi- 
bility of it being an unresolved binary cannot be ruled out. 
Thus at this stage we do not discard it as a contaminant. 

The approximate M/L transition colour is indicated 
at the bottom of the plot (estimated us ing colours from 
ILeggett et aDT l998 and Knapp et alJl2004h . and shows that 
most of our faint Praesepe candidate members should be 
very late M dwarfs. At least two however (WFC 11 and 
WFC60), have colours consistent with being early-mid L 
spectral types. These are labelled in Figure WFC60 was 
flagged as an unresolved binary in §5.2, and in the JK CMD 



it is at least 0.5 magnitudes brighter than candidate mem- 
bers with similar J-K colours. This candidate could be an 
unresolved binary L dwarf in the cluster. 



5.4 Contamination 

In this section we estimate the expected level of field con- 
tamination at the faint end of the Praesepe sequence. We 
do this for three K- band magnitude ranges: 15.5-16.0, 16.0- 
16.5 and 16.5-17.0. In the cluster member region of the IK 
CMD, it can be seen that these ranges correspond to I-K 
colours of ~ 3.5, 4.0 and 4.4 respectively. 

These colours relate to spectral types ~ M6, M8 
and L0 respectively. The absolute magnitude s of field 
stars with these sp ectral types (from iDahn et al.ll20o3 and 
IKnapp et aT] r2004) means those at 160-200 pc will overlap 
the cluster sequence in a CMD. 

The sky area covered for each magnitude range varies, 
since the different surveys considered have different photo- 
metric depth. All of the surveys considered (excepting the 
HSHJ survey) contribute to the brightest magnitude range. 
However, the Iz survey does not contribute to the two fainter 
ranges. All except the HSHJ and Iz surveys contribute to 
the middle magnitude range, and only the WFC and M98 
surveys contribute to the faintest magnitude range. After 
accounting for overlap between certain survey areas, we find 
that the brightest magnitude range covers a volume of ~ 1.2 
x 10 3 pc 3 , and the two fainter ranges each cover a volume 
of 1.0 x 10 3 pc 3 . 

Luminosity functions predict late M and early L field 
dwarf s densities of $ = ~ 2.5 x 10~ 3 and 1.9 x 10~ 3 s tars 
pc" 3 jKirkpatrick et alJ Jl994l) and ICruz et al.1 <2003ft re- 
spectively). We hence expect a field contamination of ~ 3 
late M dwarfs and ~ 2 early L dwarfs amongst the candi- 
dates. 

While clearly the levels of contamination amongst our 
candidate members with the late M spectral types is fairly 
low (~20%), it is entirely possible that both of our candi- 
date L type members could be field stars. A more rigorous 
assement of their membership status can be performed by 
measuring their proper motions. 

5.5 J-H,H-K 2 colour diagram 

Figure |S| shows our candidates in the J-H,H-K 2-colour 
diagram. The symbols are the same as in the CMDs. 
NEXTGEN (dashed lines) and DUSTY (solid lines) model 
colours are over-plotted, with the surface gravity values 
shown. The bluer synthetic H-K colours for H-K<Q.2 are 
thought to be due to inaccuracies with water vapour mod- 
elling in the H band. Approximate la te M and L spectra l 
clas s locations are also indicated f from ILeggett et al.l <ll998l) 
and IKnapp et al~ll2004l1 . 

Visual inspection of the diagram shows that, as ex- 
pected, the 2 likely background stars (WFC 22 and 30) have 
colours consistent with late M dwarfs. WFC57, probably a 
galaxy, has NIR colours similar to a late L dwarf. The pos- 
sible unresolved L dwarf binary (WFC60) and the even red- 
der single candidate member (WFC11) have NIR colours 
consistent with being early-mid L dwarfs. Similarly, Rizlll 
appears too blue for an early L dwarf, and is therefore con- 



Crossing into the substellar regime in Praesepe 7 




-i • . WFC 30 

i 

ii ' 

-a 
i 

0.2 i i i I i I I I I i i I i i I I I I i I i i i 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 

H-K 

Figure 6. The Praesepe J-H,H-K 2-colour diagram. Symbols are the same as in figure [I] Spectral type estimates from lLeggett et alj 
(2002) are overplotted. 

feature. The Dusty models do not consider grains with sizes 
greater than ~ 0.24 /ira, and consequently predict negligi- 
ble grain opacities in the NIR, and no gap. Grain formation 
models however, suggest that grains could have sizes up to 
tens of microns dCooper et alj T2003 ) , which could result in 
a significant opacity increase due to Mie scattering. 

We have measured the extent of Praesepe's "M dwarf 
gap" in the I-,J-,H- and K- bands. For each band, the gap 
size was calculated as the magnitude difference between the 
faintest candidate above the gap, and the brightest candi- 
date below, ensuring that only the candidates with full IJHK 
photometry were used. The results are tabulated (with the 
corresponding values for the Pleiades from P03) in Table 
Q For all passbands, quoted errors were derived by combin- 
ing the photometric uncertainties of the two candidates that 
were used to measure the gap sizes. 

It can be seen that the gap is smaller for Praesepe than 
for the Pleiades in all the bands considered, although there is 
a small overlap if the uncertainties are taken into account. 
This is c onsistent with the i nterpretation just mentioned, 
since (as iDobbie et alJ (|2002j) explains) one expects a shal- 
lower luminosity-Te/ / relation for older populations (such as 
Praesepe), because the radius-mass relation will be steeper 
for younger populations (like the Pleiades), where objects 
are still in the relatively early stages of contraction. The 



sidered as a probable non-member. It is shown in Figure [b] 
as an open square. 



6 THE PRAESEPE SEQUENCE 
6.1 M dwarf gap 

In FigureQwe show a zoomed in portion of the Praesepe IK 
CMD with non-members omitted. The four binary tracks 
(shown with lines joined by plus signs) will be discussed in 

JO 

It can be seen t hat there is a dearth of cluster candidates 
from if=15.3-15. 5. IDobbie et alJ J2002ft and P03 presented 
photometric evidence for a similar paucity of M7-M8 spec- 
tral types in other clusters, as well as in star forming regions 
and the field. This paucity has been referred to as the "M 
dwarf gap" (P03). The fact that the M dwarf gap has been 
observed in both young associations such as the Pleiades and 
9 Orionis as well as amongst older populations (M4 and the 
field), suggests that its origins are due to a sharp local in- 
crease in the slope of the magnitude-mass relation for late-M 
spectral types (T e // ~ 2700-ftT). At this T e //, models predict 
that dust g rains begin to con dense in t he outer atmo s pheric 
layers (e.g. iTsuii et alJll996l) . and as IDobbie et~aH (bo02fi 
suggests, grain opacities could be responsible for the gap 
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Table 1. Sizes of the "M dwarf gap" in the I-,J-,H- and K- bands 
for the Pleiades and Praesepe 



Cluster 


I 


J 


H 


K 


Pleiades 


0.5±0.05 


0.3±0.06 


0.3±0.04 


0.3±0.05 


Praesepe 


0.30±0.03 


0.25±0.05 


0.24±0.05 


0.23±0.03 



relative extent of the Praesepe M dwarf gap in the four pho- 
tometric bands considered is comparable for both clusters. 
We will, however, discuss more subtle colour changes across 
the gap in the next section. 



6.2 The faint single star sequence 

In order to study how the bulk properties of the Prae- 
sepe candidates change for K >15, we split the Praesepe 
candidate members into the following A"-band magnitude 
bins: bin A (15.0 - 15.3) represents the candidates slightly 
brighter than the M dwarf gap (discussed in the last section) . 
Bins B and C (15.5-15.75 and 15.75-16.00 respectively), are 
below the gap. Bin D covers all the fainter candidates (ex- 
cept WFC 11) down to K =16.8. These 4 bins are indicated 
by dashed lines in Figure [7| 

In order to determine an "observational photometric se- 
quence" , we used the NEXTGEN model isochrone above the 
gap region in the CMD (since this model agrees well with 
the data in this range). To bridge the M dwarf gap we have 
used our measured gap sizes (see Table 1) to provide colour 
and magnitude changes from K=15.3 - 15.5. Below the gap 
(bins B-D), we have averaged the colours and magnitudes 
of our single star candidate members in each bin. Our "ob- 
served cluster sequence" is shown as a solid line in Figure 

HI 

When comparing the observed sequence to the 
model isochrones, it is apparent that for I -K =3.3-4 the 
NEXTGEN model is still quite reasonable. However, for 
I-K><i the candidates appear to drop slightly below 
the NEXTGEN is ochrone, moving clos er to the DUSTY 
isochrone. P03 and I Jameson et all <l2002f) pointed out a sim- 
ilar transition from the NEXTGEN to DUSTY isochrones 
for the Pleiades cluster sequence from 7-A'=4.1-4.5. 

We trace the colour changes of the single source can- 
didates (ie. avoiding sources flagged as unresolved binaries) 
through the 4 K-band magnitude bins by indicating where 
the bin members lie (with boxes) on a zoomed in version of 
the 2-colour diagram (shown in Figure [HJ. The boxes were 
constructed such that single candidate members from each 
of the if-band magnitude bins are all contained within the 
respective box dimensions in the diagram. 

Bin A, which is populated by candidates above the M 
dwarf gap, appears to be slightly bluer in (J-H) than the 
NEXTGEN model colours, lying closer to the DUSTY model 
colours. We expect no dust condensation at these T e //s, so 
this difference probably results fr om incompletenes s in the 
NEXTGEN H 2 opacities (eg. see I Jones etaDl200^ . Bin B 
is slightly redder in H-K than Bin A, but has similar J-H 
colours. This colour change could be caused by the addition 
of new grain opacities, consistent with the appearance of at- 
mospheric dust across the M dwarf gap. However, note that 



no colour changes were observed across the M dwarf gap 
in the Pleiades (see P03). If the M dwarf gap does result 
from the first appearance of atmospheric dust grains, then 
they must have different effects for different surface gravities. 
This would not be unexpected since we would expect grain 
growth and rain-out to be affected by both surface grav- 
ity and pressure. The fainter bins (C and D) have colours 
consistent with the Dusty models. 



7 THE MASSES OF LOW MASS CANDIDATES 

We have estimated the masses of all single Praesepe candi- 
dates that appear in our candidate sequence region (see Fig- 
ure 4). The exception to this is hill, which was identified 
in H5.5I as a likely non-member. Masses were estimated us- 
ing the NEXTGEN and Dusty models (with (m-M) =6.16) 
for isochrone ages o f 0.5 and 1 Gyrs, and solar metallicity. 
iMartfn et alJ ((2000:) note that mass estimates from J- band 
magnitudes are the least sensitive to the choice of model, ei- 
ther NEXTGEN or DUSTY. Therefore, we choose to derive 
masses from this passband where possible. Where J- band 
photometry was absent, masses were estimated from the K- 
band magnitudes. 

Linear interpolation was used in between the isochrone 
mass points. The DUSTY mass estimates are slightly higher 
than the NEXTGEN ones, and, rather more obviously, the 
younger mass estimates are lower than the older ones. The 
mass estimates are given in Table lA"3l Estimates taken from 
the K- band are slightly higher than if taken from the J- 
band. 

It can be seen from the table that we have identified at 
least 15 single VLMS candidates (with masses ^ 0.075 Mq). 
The five faintest candidates could be sub-stellar or stellar 
depending on the age or model we adopt. For example, four 
of these are substellar if a 0.5 Gyr cluster age is assumed, 
regardless of which model is used. 



8 PRAESEPE BINARIES 
8.1 Binary mass ratios 

In Figure [7| we have over-plotted 6 binary tracks, which we 
use to estimate mass ratios of our Praesepe binaries. We 
calculated these tracks by combining single star points from 
the cluster sequence. Each binary track begins at some single 
star point, and the affect of adding an unresolved companion 
is simulated by combining the photometric brightness of this 
single star point with other single star points from lower 
mass points in the sequence (beginning with the lowest mass 
point on the sequence, and then increasing the companion 
mass until an equal mass binary is reached, 0.75 magnitudes 
above the original single star point) . 

When estimating binary mass ratios (q), one must ac- 
count for the cluster depth which can make sources appear 
brighter without the need for a binary companion, as well 
as phot ometric uncertainty. The tidal radius of Praesepe is 
12.1 pc llHolland et al]l200(fl giving a depth effect ~ ± 0.15 
magnitudes. Photometric uncertainty is typically ±0.2 in 
I-K. Mass ratios were estimated by comparing unresolved 



Crossing into the substellar regime in Praesepe 9 




Figure 7. K, I-K CMD. Symbols are the same as for figure^ Solid line shows the empirical sequence that tracks the move away from 
the NEXTGEN isochrone and onto the dusty sequence. 0.5 Gyr NEXTGEN and Dusty models are shown as dotted and dashed lines 
respectively. 4 binary tracks are also shown. 





Figure 8. J-H,H-K 2-colour diagram. Symbols are the same as for figure [I] Dotted lines define the perimeters of 4 K magnitude bins, 
discussed in the text. 
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Table 2. Binary mass ratios assuming a cluster depth of ± 0.15 r ' 
using 0.5 Gyr and 1.0 Gyr NEXTGEN and Dusty models. 



Table 3. Praesepe binary fractions. Mass and q ranges estimated 
using Tables 2 and A3. 





q (0.5 Gyr) 


q (1.0 Gyr) 




NG 


DUSTY 


NG 


DUSTY 


RIZ2 


0.8-1.0 


0.8-1.0 


0.9-1.0 


0.9-1.0 


RIZ21 


0.6-1.0 


0.5-1.0 


0.6-1.0 


0.6-1.0 


RIZ18 


0.7-1.0 


0.6-1.0 


0.6-1.0 


0.6-1.0 


IZ126 


0.4-1.0 


0.4-1.0 


0.5-1.0 


0.5-1.0 


WFC169 


0.7-1.0 


0.7-1.0 


0.8-1.0 


0.7-1.0 


WFC118 


0.6-1.0 


0.6-1.0 


0.7-1.0 


0.6-1.0 


WFC88 


0.7-1.0 


0.7-1.0 


0.7-1.0 


0.6-1.0 


WFC60 


0.8-1.0 


0.7-1.0 


0.8-1.0 


0.8-1.0 



binary candidates to the binary tracks in the CMD, and al- 
lowing for both these forms of uncertainty. The results are 
tabulated in Table [5] 

The q values range from 0.4-1.0. This represents a 
rather larger range of q values than was found for Pleiades 
BDs by P03 (q=0. 7-1.0). However, the photometric uncer- 
tainties are slightly larger for the Praesepe candidates, and 
the very low mass unresolved binary populations identified 
in the two clusters cannot be shown to be significantly dif- 
ferent (in their q range) by this analysis. 

8.2 The binary fraction 

P03 determined binary fractions (BFs; defined as the num- 
ber of binary systems divided by the total number of sys- 
tems in some photometric range) for Praesepe in 4 colour 
bins for 1 ^ I-K ^ 3.6, corresponding to an approximate 
mass range 1.0 - 0.09 Mq. We extend this to lower masses 
by considering two redder bins from I-K =3.3-3.9 and 3.9- 
4.6. For the first bin the Iz survey is not complete over the 
full single star sequence. Therefore, to avoid underestimat- 
ing the relative number of single star sources we only con- 
sidered Riz and WFC candidates. With 5 unresolved binary 
and 11 single source candidates, this bin has a BF~31%. 
For the second bin, the Riz survey is not complete over the 
full single star sequence. However, the Riz candidates are 
all contained within the WFC survey area, and would have 
been discovered by this survey had they not already been 
identified. We therefore considered all sources in this bin, 
resulting in 3 unresolved binary and 7 single source candi- 
dates, or a BF~30%. These results are summarised in Table 
[3] where uncerta inties have been calcu lated assuming bi- 
nomial statistics l)Bureasser et al.ll2003h . Note that the low 
number statistics give large uncertainties, and that the BF 
for the redder bin is a lot more uncertain even than this, 
since several of these candidates could be non-members (see 
Sections 5.4 and 10.4). 

The Praesepe BFs are displayed in Figure [5] with the 4 
higher mass B Fs from P03 o ver- plotted. Also shown are the 
BFs from lClose et alJ J2003I) and lBurgasser et alJ (120031) for 
VLMS and BD binaries in the field with separations >1AU. 
Although our lowest mass point is too uncertain to accu- 
rately constrain the BF, the higher mass point shows that 
the BF is decreasing as one goes below O.IMq. Even if our 
expected level of contamination (~3 late M dwarfs) removes 
only single star candidates (thus pushing up the BF) , the BF 
would only go up to 391^%, which would still suggest a BF 



I-K range 


MASS (0.001 M ) 
NEXTGEN DUSTY 


q 


13 F 




0.5 


1.0 


0.5 


1.0 






3.3-3.9 
3.9-4.6 


89±9 
73±7 


90±5 
81±5 


70±5 
63±4 


78±4 
74±5 


0.5-1.0 
0.4-1.0 


31±J a 
30±« 



Figure 9. Praesepe binary fractions. Two points from the liter- 
ature are displayed as open diamonds. 



decrease. This decreasing BF trend appears to be consis- 
tent with the lower BFs seen for field sources. P03 however, 
observed an increasing unresolved BF into the Pleiades BD 
regime, reaching 50ii° f° r O.O7M0. This suggests that dif- 
ferent clusters could have different very low mass BFs. Such 
differences could result from different cluster environments, 
and in this context the distribution of different cluster BFs 
could provide an important observational test of VLMS and 
BD formation theories. However, before any firm conclusions 
are drawn, it is clearly important to improve the robustness 
of the measured cluster BFs by confirming (or not) cluster 
membership from proper motions. 



9 LUMINOSITY AND MASS FUNCTION 

We divided our cluster members into 3 magnitude ranges 
covering K= 15-15.8, 15.8-16.3 and 16.3-17. For each mag- 
nitude range, cluster candidates were counted in 0.5 degree 
rings from the centre, out to a radius of 2.5 degrees for the 
brightest bin, and out to 2.0 degrees for the 2 fainter bins. 
We did not consider candida tes found in th e region of the 
sub-cluster identified bv lHolland et al] ll200d) . 3 pc from the 
centre of Praesepe. This sub-cluster could be a smaller older 
open cluster that has collided with Praesepe, and we avoided 
candidates within 1 degree of the sub-cluster centre. The 
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positions of candidates in the magnitude bins A"=15-15.8, 
15.8-16.3 and 16.3-17.0 are shown in figure 1101 as crosses, 
squares and asterisks respectively. It can be seen from the 
figure that there is an over-density of faint objects in the 
sub-cluster region. It is clearly important to isolate (and ig- 
nore) this region when measuring the LF and MFs of the 
main cluster, since the s ub-cluster is though t to be older 
than the rest of Praesepe jHolland et all2000T) . and it would 
otherwise be possible for faint low-mass sub-cluster stars to 
masquerade as Praesepe BDs. 

Ring counts were then converted into surface densities 
by dividing by the appropriate survey areas within each ring 
(ignoring the sub-cluster region) , and expected levels of field 
star contamination (see Section were subtracted off the 
surface density profiles assuming a uniform spatial distribu- 
tion over the cluster. We then assumed that the cluster could 
be re presented by a King surface densit y distribution ( King 
1962) (as has been don e previously bv|Pinfield et al]|200Ct 
iJameson et"alll2002l and lHolland et al.ll200fj) . This function 
takes the form: 



fs = k 



y/l + (r/r c y y/l+(rt/r c ¥ 



(1) 



Here, f s is the surface density, r the radius from the clus- 
ter centre, k a normalization constant and r c and rt the 
core and tidal radii respectively. Assuming r t = 12.1 pc 
jHolland et alJ |2000) , and using Poisson uncertainties for 
the surface density profiles, we obtained the best fit value 
of r c using the whole sample, and found r c =7.56 pc. This 
was then assumed to be constant for each magnitude bin, 
so we proceeded to minimize the x 2 statistic of the King 
profile to the 3 observed density distributions. We then de- 
termined the total number of cluster stars in each bin using 
the surface integral (out to rt) of equation 



irr c K < ln(l + x t ) 



■ x t 



1)(VT 



■ Xt 



1+Xt 

(2) 

where Xt = (r t /r c ) 2 . Our results are summarised in Table 
2] The resulting luminosity function ($x) is shown in figure 
llll as filled diamonds. Also shown is the luminosity function 
from HSHJa (converted into &k), although note that the 
faintest three HSHJ points (open circles) are lower limits 
due to photometric incompleteness. 

The mass function then follows directly from the lumi- 
nosity function: 



N(m) = $j 



dK 
dm 



(3) 



We estimated 4^- for our 3 magnitude ranges using both the 
0.5 Gyr and l.o"byr NEXTGEN and DUSTY isochrones to 
define the mass limits of each bin, and derived the cluster 
mass function per O.IMq interval (see Table0]l. 

The log-norma l Pleiades mass function as derived by 
lAdams et al.l (|20Jll|) is displayed on Figure 1121 as a dashed 
line, and is normalised to the Praesepe function at 0.5 
M(7). An upper lim it Hyades mass function point from 
iDobbie et all 120021 in the range 0.08-0.06 M is also shown 
for comparison. 

It can be seen from Figure 11 that both the Praesepe 
mass function derived here and the Hyades mass function 




Figure 10. Sky plot showing the positions of single candidates 
fainter than K-=15.0. Crosses, squares and asterisks represent 
single candidates in the magnitude bins K=15— 15.8, 15.8-16.3 
and 16.3-17.0 respectively. Annuli of radius 0.5°, 1.0°, 1.5° and 
2.0° from the centre of the cluster are also shown as dotted lines. 
The position of the excluded subcluster is represented by the solid 
line. 




Figure 11. Praesepe iC-band luminosity function. Diamonds are 
points derived in this work. Circles are from the HSHJ survey, 
where the last three HSHJ points are lower limits. 
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Table 4. Luminosity and mass function bins with derived values as discussed in the text. 















Mass Ranf 


;e (0.001 M Q ) 


N(m) /0.1 M 


K Range 


Surveys 


N°- 


Contam. 


N(r t ) 




0.5 Gyr 


1.0 Gyr 


0.5 Gyr 


1.0 Gyr 


15.0-15.8 


Riz,Iz,M98,WFC 


15 


3.0 


34.6±8.9 


43.2±11.2 


87-123 


91-122 


96.0±24.8 


111.5±28.8 


15.8-16.3 


Riz,M98,WFC 


3 


2.5 


3.3±1.9 


6.6±3.8 


74-93 


84-96 


17.3±10.0 


27.4±15.8 


16.3-17.0 


M98,WFC 


3 


2.0 


8.8±5.1 


12.6±7.3 


61-78 


74-86 


51.9±29.9 


73.5±42.4 



1000 i i r 




1.0 0.1 
Moss 

Figure 12. Praesepe mass function. Symbols ar e same as for 
figure [TT1 The log-normal Pleiades MF best- fit from lAdams et all 
(2001) is shown as a dotted line, normalized at 0.5 Mq. Also 
shown is a 0.06-0.08 Mq MF u pper limit for the Hyades cluster 
taken from lDobbie et al.ll20o3) . 

falls below that of the Pleiades at ~0.1M e . The Hyades and 
Praesepe are thought to be of a similar age (~ 0. 6 Gyr), com- 
pared to the younger 0.13 G yr Pleiades cluster <|Burke et alJ 
1200^ . This age gap has lead lDobbie et alJ J2002ft to suggest 
that the paucity of low-mass Hyads is caused by dynamical 
evolution of the environment. In this scenario, the lowest- 
mass members of an association are, over time, preferentially 
ejected from th e cluster - a theory supported by the jV-body 
simula tions of Ide la Fuente Marcos fc de la Fuente Marcos! 
(2000) in which a 0.6 Gyr open cluster has retained at best 
only ~ 17 percent of the BDs that were originally bound to 
the association. Evidently, our Praesepe mass function is in 
strong agreement with this scenario. 

Despite the qualitative similarities, however, between 
the mass functions of Praesepe and Hyades, there are to 
date no Hyad substellar candidates, in comparison to the 
Praesepe candidates proposed in this work. As a result, the 
Hyades low-mass MF upper limit is only marginally consis- 
tent with our lowest mass Praesepe point, which suggests 
that the Praesepe MF may be somewhat higher than that 
of the Hyades in this mass range. Clearly, differing dynami- 
cal processes would not be an issue here since both clusters 
are of a similar age and can be assumed therefore to be at a 
similar stage of evolution. If this is a real feature of the MF, 



rather than as a consequence of low-number statistics, we 
speculate that this difference can perhaps be explained by a 
substellar mass function sensitivity to natal molecular cloud 
properties such as densi ty, metallici t y, tem perature etc. In- 
deed, the simulations of lElmegreenl (|200(J) suggest that the 
shape of the MF below the observed 0.1 Mq turnover is 
dependent upon the physical cloud properties, in contrast 
to the universal Salpeter IMF that characterises the inter- 
mediate to high mass domain. We stress, though, that it is 
imperative to increase the known substellar population in 
a variety of associations before this hypothesis can be fully 
tested, and that the possibility of a common substellar MF 
in both environments cannot be ruled out. 



10 ASTROMETRY 

We have mentioned, in previous sections, how important 
proper motions are in confirming cluster membership. The 
peculiar motion of Praesepe is large compared to that of 
typical field stars, and the intrinsic velocity dispersion of the 
cluster is small (HSHJb) . With a base-line of a few years, it 
is thus possible to measure the proper motion of candidates 
with sufficient accuracy to separate true cluster members 
from field star contamination. In this section we describe 
how we have measured proper motions for a subset of the 
Praesepe candidates we consider in this work. 

10.1 Epoch data 

First epoch data for the astrometric analysis of Praesepe RIZ 
candidates was taken in 1993 using the 2.5 m Isaac Newton 
Telescope (INT), providing a baseline of ~ 8.9 years. For 
Roque PR1, we use the original image acquired in Febru- 
ary 1996, also taken with the INT, with a baseline of ~ 5.9 
years. We use our WFC images as 2nd epoch data. The im- 
ages in which RIZ 21 and 23 were recovered were acquired 
on the third night of the observing run, which could not be 
photometrically calibrated due to poor conditions. Conse- 
quently, these two candidates are not plotted in Figure |3] 
However, we do not expect this to affect the astrometric 
measurements. For all images, we use only the / band for 
the astrometry. 

10.2 Measurements 

The process of deriving the astrometry of the cluster can- 
didates is as follows. The pixel positions of the candidates 
on both epoch images (x lst ,yist),{x2nd,y2 n d) are first accu- 
rately determined using the IRAF: CENTRE routine. The 
positions of presumed background point sources were also 
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measured in both epochs, to be used as astrometric refer- 
ences. IRAF:XYXYMATCH was used to calculate the trans- 
formation function which maps the two epoch images. GE- 
OMAP and GEOXYTRAN were then used to transpose the 
second epoch coordinates of the candidates into the first 
epoch coordinate system. Any reference sources lying more 
than a Gaussian 3 a from the median of absolute deviation 
were rejected and the transform redetermined. 

The relative proper motions along the x and y axes p x , y 
of the CCD then follows from the pixel displacement: 



px 



(AX) 
At 



At 



(4) 



PY = 



(AY) _ Y 2 , 



(5) 

were found to 



At At 

The uncertainties of the transformations a x , y 
be « 0.2 pixels, corresponding to proper motion errors of « 
10-15 masyr -1 , in both RA and Dec direction. 

10.3 Membership probabilities 

The metho d of estab l ishing membership probabilities is that 
outlined bv lSandersI il97lT) . and utilizes the following likeli- 
hood function: 



(6) 



where 



2ttH x E v (1- P 2 



(7) 

here, n c is the normalized number of cluster stars, 
and p is a correla tion coefficient. Following the method of 
IWang et alJ <ll995f) , a follows from the proper motion uncer- 
tainties a x , y and an intrinsic variance ctq : 



2 2 2 

O = + a x,y 



(8) 



The distribution functions of the cluster and field & c , <!>/ 
are evaluated as: 



$c(Hx,p, y ) = exp 



and 



& f (p x ,p y ) = exp 



Px pxc \ 2 . / py l^yc ^ 

V a 



Px ~ Pxf 



(9) 



2(1 -p") 



2p{p x - Pyf)(Py ~ Pyf) ( Py - Pyf 



(10) 



here, a is the proper motion dispersion, p xc ,p yc are the 
cluster mean proper motions, and p x f,P y f are the field mean 
proper motions with standard deviation s Ti x and T, v . 

The free parameters were found by[Wang et alJ (Il995h 
using a conjugate gradient method using a large Praesepe 
census, and are reproduced in Tabled] 



Table 5. Praesepe maximum likelihood parameters 



parameter 


value 


pxc 


-32.88 ± 0.09 mas yr" 1 


Pyc 


-15.87 ± 0.09 mas yr" 1 


Pxf 


0.23 ± 0.08 mas yr" 1 


Vyf 


-7.18 ± 0.26 mas yr" 1 


CO 


1.45 ± 0.09 mas yr" 1 


p 


0.002 ± 0.016 




■A. 



-20 
Wos/yr 



Figure 13. Proper motion vector point diagram for objects lo- 
cated in Praesepe fields. Asterisks arc the sources used as astro- 
metric references. Filled stars show the RIZ candidates, and the 
triangle is RPrl. Field and c luster members (P > 70 percent) as- 
signed in lWang et alJ Il995f) are shown as small and large filled 
circles respectively. 



10.4 Astrometric results 

Proper motions and membership probabilities of RIZ ob- 
jects from P97 that have JHK colours consistent with clus- 
ter memberships are shown in Table El and illustrated in a 
vector point diagram (VPD) in Figure H3l as filled stars with 
error bars. We were unable to measure the proper motion of 
RIZ 4,11,18 and 24 due to insufficient shared spatial cover- 
age between epochs, which would have resulted in unreliable 
transformations. 

Proper motions from lWang et all il995h are overplotted 
as small circles, or large circles where membership probabili- 
ties were evaluated to be greater than 70 percent. These help 
to clarify the cluster position at p x , y ~ -33,-16 masyr -1 . All 
stars that were used as references for the astrometric trans- 
formation are also shown on the VPD as asterisks. As ex- 
pected, they cluster around the origin of the diagram, since 
they were selected on the basis of having a negligible motion. 

All of the RIZ candidates have an evaluated member- 
ship probability greater than 70 percent, and are therefore 
consistant with cluster membership within the astrometric 
uncertainties. Our astrometric analysis of Roque Prl, shown 
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Table 6. Astrometric results. J is also WFC91. 



Object fjt x mas yr 1 fi y mas yr 1 Prob mem 



RlZlt 


-34.45 ± 18.48 


-0.95 ± 16.14 


86.72 


Y 


RIZ2 


-41.87 ± 16.37 


-9.44 ± 12.23 


94.33 


Y 


RIZ8 


-31.80 ± 12.92 


-11.48± 9.37 


97.87 


Y 


RIZ20 


-27.80 ± 16.52 


-7.64 ± 15.66 


86.56 


Y 


RIZ21 


-38.75 ± 16.37 


-1.56 ± 13.05 


72.11 


Y 


RIZ23 


-47.12 ± 19.75 


-2.83 ± 12.77 


75.41 


Y 


RPrl 


-33.55 ± 10.13 


13.67 ± 13.52 


51.24 


N? 



as a triangle on the VPD, suggests that this object is most 
likely a non-cluster member. RPrl is more likely to be a field 
star in the line of sight to the cluster. 



11 CONCLUSIONS & FUTURE WORK 

In this paper, we have presented the results of a deep optical 
survey (I~21.3) and NIR follow-up of 2.6 square degrees of 
the Praesepe open cluster. We have used our census of cool 
red members to place constraints on the bottom end of the 
cluster mass function. Our results indicate that the dynam- 
ical evolution of the cluster is likely leading to the evapo- 
ration of the lowest mass members. However, the relative 
depletion of the very-low mass stellar and substellar mem- 
bers may be less than that reported for the similarly aged 
Hyades cluster, suggesting that perhaps the cluster MF is 
sensitive to natal molecular cloud properties. 

To test the robustness of these arguments, it is imper- 
ative to ascertain beyond doubt the validity of the candi- 
dates' claims to cluster membership. We have done this us- 
ing proper motion measurements for several low-mass Prae- 
sepe candidates, but our priority is to continue this for all. 
Only then can we probe with confidence the substellar MF. 
Furthermore, a bona- fide census would strengthen the com- 
parison between low-mass binary fractions in this cluster 
compared to earlier clusters such as the Pleiades. Our re- 
sults suggest a difference between the two environments, but 
its cause or quantitative nature cannot be elaborated upon 
without further census expansion. 

The photometric sample we have compiled have re- 
vealed a paucity of candidates close to the region where the 
main sequence begins to deviate away from the NEXTGEN 
track. This feature has been observed in younger regions 
too, and would appear to originate because of a steepening 
in the mass-luminosity relation at late M spectral type. We 
propose that further measurements, both photometric (e.g 
L- band) and spectroscopic, be made of candidates which 
are located just above and below the gap in both this clus- 
ter and others to isolate the underlying physics, and to test 
the sensitivity of dust formation in the atmosphere of early 
cluster L dwarfs to the fundamental properties of their as- 
sociation. 

Finally, to further the results presented in this work 
it is of utmost importance to increase the known number 
of VLMS and BDs in young SFRs and open clusters. In 
the near future, this need will primarily be met by the ad- 
vent of automated deep surveys conducted over wide areas. 
One such effort is UKIDSS, the collective name for a se- 
ries of large-scale IR surveys using the WFCAM instrument 



on the UK Infrared Telescope. The Galactic Cluster Survey 
(GCS) in particular will survey 10 associations over the next 
few years on an unprecedented scale, uncovering many more 
cluster BDs that will go a long way to addressing some of 
the issues raised in this work. 
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APPENDIX A: TABLES OF PHOTOMETRY, 
MEMBERSHIP AND BINARITY CRITERIA, 
AND MASS ESTIMATES 
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Table Al. Photometry of the WFC candidates that were initially considered interesting (see Section 3). Final photometric membership 
criteria are given in the last 3 columns. J and f are also Rizl and 2MASSJ08350622+1953050 respectively. 



Membership criteria 

WFC RA (J2000.0) Dec. (J2000.0) I c J H K IK JK JHK 



11 


8 45 


41.05 


19 


38 


02.6 


21.37 


18.37 ±0.04 


17.51 ±0.02 


16.84 ±0.02 


Y 


Y 


Y 


22 


8 44 


34.29 


19 


14 


23.3 


21.20 


18.37 ±0.04 


17.74 ±0.06 


17.22 ±0.04 


N? 


N? 


? 


24 


8 38 


32.13 


19 


59 


27.1 


21.19 


17.99 ±0.03 


17.30 ±0.02 


16.78 ±0.02 


Y 


Y 


Y 


30 


8 41 


56.61 


19 


20 


34.8 


21.13 


17.91 ±0.04 


17.45 ±0.03 


17.03 ±0.04 


N? 


N 


? 


53 


8 46 


39.19 


19 


27 


37.0 


20.90 


17.83 ±0.03 


17.18 ±0.03 


16.64 ±0.03 


Y 


Y 


Y 


57 


8 37 


19.80 


19 


38 


51.4 


20.79 


18.82 ±0.06 


17.90 ±0.06 


16.91 ±0.03 


N? 


N 


N 


60 


8 41 


25.03 


19 


09 


25.0 


20.79 


17.61 ±0.02 


16.80 ±0.02 


16.22 ±0.01 


Y 


Y 


Y 


76 


8 37 


27.61 


19 


37 


33.0 


20.51 


17.58 ±0.01 


16.83 ±0.01 


16.46 ±0.02 


Y 


Y 


? 


81 


8 41 


50.21 


19 


06 


18.4 


20.33 


17.47 ±0.01 


16.80 ±0.01 


16.30 ±0.01 


Y 


Y 


Y 


87 


8 38 


02.34 


19 


46 


09.0 


19.83 


17.14 ±0.01 


16.64 ±0.01 


16.06 ±0.01 


Y 


Y 


? 


88 


8 40 


43.86 


19 


06 


00.8 


19.83 


_ 


_ 


15.59±0.02 


Y 


_ 


_ 


91* 


8 36 


54.46 


19 


54 


15.3 


19.75 


17.07 ±0.01 


16.47 ±0.01 


15.99 ±0.01 


Y 


Y 


Y 


94 


8 35 


51.93 


19 


20 


33.6 


19.69 


17.04 ±0.01 


16.42 ±0.01 


15.98 ±0.01 


Y 


Y 


Y 


103 1 " 


8 35 


06.12 


19 


53 


05.3 


19.56 


_ 


_ 


14.28±0.01 


N 


_ 


_ 


104 


8 41 


47.75 


19 


48 


37.7 


19.56 


_ 


_ 


15.80±0.02 


Y 


_ 


_ 


115 


8 36 


44.74 


19 


17 


23.9 


19.38 


_ 


_ 


15.93±0.03 


Y 


_ 


_ 


118 


8 37 


38.73 


19 


03 


09.6 


19.30 


_ 


_ 


15.53±0.02 


Y 


_ 


_ 


135 


8 38 


07.96 


19 


43 


28.6 


19.11 


_ 


_ 


15.74±0.02 


Y 


_ 


_ 


136 


8 45 


39.10 


19 


38 


04.2 


19.09 


- 


- 


15.80±0.02 


Y 


- 


- 


138 


8 36 


52.03 


19 


00 


42.2 


19.08 






15.78±0.08 


Y 






140 


8 36 


21.20 


19 


38 


45.5 


19.06 






15.56±0.02 


Y 






142 


8 34 


54.28 


19 


02 


32.7 


19.03 






15.62±0.02 


Y 






150 


8 36 


48.06 


19 


49 


02.2 


18.92 






15.52±0.02 


Y 






159 


8 37 


41.90 


19 


09 


31.0 


18.82 






15.59±0.02 


Y 






163 


8 40 


53.59 


19 


40 


58.3 


18.78 






15.60±0.02 


Y 






169 


8 36 


55.63 


19 


36 


15.5 


18.69 






14.99±0.01 


Y 






191 


8 38 


11.96 


19 


59 


43.5 


18.25 






15.22±0.02 


Y 






193 


8 41 


27.76 


19 


07 


19.7 


18.21 






15.27±0.02 


Y 






200 


8 46 


16.82 


19 


52 


30.2 


18.17 






15.05±0.02 


Y 






204 


8 36 


29.90 


19 


38 


37.5 


18.08 






15.14±0.02 


Y 






207 


8 39 


39.37 


19 


47 


54.9 


18.04 






15.19±0.02 


Y 






208 


8 40 


50.31 


18 


23 


59.8 


18.02 






14.82±0.01 


Y 






210 


8 38 


50.24 


19 


30 


18.6 


18.02 






15.14±0.02 


Y 






216 


8 46 


48.87 


19 


32 


09.4 


17.93 






14.98±0.02 


Y 






220 


8 37 


02.23 


19 


52 


07.6 


17.88 






14.76±0.01 


Y 
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Table A2. HSHJ photometric non-members. HSHJ{2,8,9,11, 39,54,73,88,109, 110,217,413} had no 2MASS counterpart « 3"). 
HSHJ{59,114,294} were flagged in 2MASS as contaminated, due most likely to either photometric confusion from nearby sources or 
to diffraction spikes. HSHJ{509,510,515} have red colours typical for foreground field dwarfs but, in light of the fact that they all appear 
on the same HSHJ plate, we consider their red colours as potentially suspect, and do not flag them as non-members. 



HSHJ RA (J2000.0) Dec. (J2000.0) I c K 



46 


8 33 


43.55 


19 


56 


34.94 


17. 


26 


14.94 


±0.10 


98 


8 36 


05.97 


22 


26 


15.98 


12. 


62 


11.89 


±0.02 


174 


8 37 


35.96 


18 


41 


25.25 


16. 


78 


14.55 


±0.07 


199 


8 38 


16.04 


21 


23 


29.87 


12. 


63 


11.44 


±0.02 


200 


8 38 


14.54 


19 


52 


34.41 


15. 


17 


13.17 


±0.03 


208 


8 38 


36.26 


22 


32 


27.56 


12. 


62 


11.64 


±0.02 


231 


8 38 


52.07 


18 


07 


44.53 


12. 


89 


11.61 


±0.02 


315 


8 40 


21.35 


19 


10 


53.35 


12. 


20 


10.75 


±0.03 


323 


8 40 


33.53 


19 


38 


00.48 


11. 


16 


10.17 


±0.02 


332 


8 40 


53.87 


19 


59 


55.45 


13. 


99 


12.10 


±0.02 


336 


8 40 


55.69 


19 


51 


51.84 


11. 


.88 


10.79 


±0.02 


387 


8 41 


44.13 


22 


04 


17.41 


16. 


11 


13.98 


±0.06 


391 


8 41 


47.76 


19 


24 


43.72 


11. 


00 


10.09 


±0.02 


437 


8 43 


04.29 


19 


17 


33.76 


11. 


62 


10.56 


±0.02 


469 


8 43 


51.30 


19 


23 


54.15 


17. 


12 


14.80 


±0.12 


477 


8 44 


17.16 


19 


51 


55.64 


12. 


77 


11.12 


±0.02 


480 


8 44 


18.25 


18 


43 


50.74 


12. 


71 


11.48 


±0.02 


487 


8 44 


29.75 


19 


28 


03.43 


11. 


40 


10.59 


±0.02 



Table A3. Praesepe photometric candidates found in our cluster sequence region (see text) in a scending I-K. The Table includes WFC 
candidates (from this work), Riz and Iz candidates from P03, and RPrl iMagazzu et aljll998l) . Mass estimates (0.001 Mq) are given 
using both 0.5 and 1.0 Gyr relations, with NEXTGEN and Dusty models providing the lower and upper limits respectively. Where 
possible, mass estimates were made using the more reliable J- band predictions, otherwise the K- band was used. Unresolved binaries 
are flagged (see i|5.2l . 



Candidate I c K I c -K J-H H-K Mass (0.5 Gyr) Mass (1.0 Gyr) Bin 



WFC138 


19. 


.08 


15. 


.78 


3. 


.30 






93- 


-88 


97- 


-92 




WFC135 


19. 


.11 


15. 


.74 


3 


.37 






95 


89 


97- 


93 




WFC150 


18. 


.92 


15. 


.52 


3 


.40 






101 


96 


103 


-98 




WFC142 


19. 


.03 


15. 


.62 


3 


.42 






99 


92 


101 


96 




PJZ2 


18. 


.19 


14. 


.77 


3 


,42 


0.56 


0.34 










Y 


WFC115 


19. 


.38 


15. 


.93 


3. 


.44 






89- 


-83 


93 


-89 




IZ93 


19. 


.20 


15. 


.75 


3. 


.45 


0.55 


0.46 


94 


-88 


97 


92 




IZ126 


18. 


.70 


15. 


.22 


3 


.48 


0.37 


0.35 










Y 


WFC140 


19. 


.06 


15. 


.56 


3 


.50 






100 


-94 


102 


97 




RIZ23 


19. 


.06 


15. 


.53 


3 


.53 


0.59 


0.51 










Y 


IZ23 


19. 


.10 


15. 


.57 


3 


.53 


0.43 


0.47 


94 


93 


97 


96 




RIZ 11 


19. 


.47 


15. 


.84 


3 


.63 


0.43 


0.39 


88- 


87 


92 


91 




IZ36 


19. 


.31 


15. 


.66 


3 


.65 


0.51 


0.45 


90 


89 


93- 


92 




WFC169 


18. 


.69 


14. 


.99 


3 


70 














Y 
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